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Abstract 

We consider the possibility that the Higgs boson can act as a link to a hidden sector 
in the context of pure-glue hidden valley models. In these models the standard model 
is weakly coupled, through loops of heavy messengers fields, to a hidden sector whose 
low energy dynamics is described by a pure- Yang-Mills theory. Such a hidden sector 
contains several metastable hidden glueballs. In this work we shall extend earlier results 
on hidden valleys to include couplings of the messengers to the standard model Higgs 
sector. The effective interactions at one-loop couple the hidden gluons to the standard 
model particles through the Higgs sector. These couplings in turn induce hidden glueball 
decays to fermion pairs, or cascade decays with multiple Higgs emission. The presence 
of effective operators of different mass dimensions, often competing with each other, 
together with a great diversity of states, leads to a great variability in the lifetimes and 
decay modes of the hidden glueballs. We find that most of the operators considered in 
this paper are not heavily constrained by precision electroweak physics, therefore leaving 
plenty of room in the parameter space to be explored by the future experiments at the 
LHC. 
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1 Introduction 



In view of the construction of the Large Hadron Collider (LHC) and its detectors, there have 
been major efforts to develop ideas for what physics may lie beyond the standard model (SM). 
This work is often motivated by a number of problems that still remain unsolved such as the 
hierarchy problem, baryogenesis and the origin of dark matter, among others. Nevertheless, 
it is also possible that new physics is unrelated to the above phenomena and may show up in 
a way that is not predicted by existing theories. To ensure that such a physics does not go 
undiscovered requires precise understanding of how new physics will reveal itself in the next 
particle-physics experiments. Therefore, it is prudent to explore other approaches to studying 
new theories, with an emphasis in models with varied experimental signatures. 

The hidden valley scenario [U El EJ HI G3 E] has been proposed as an interesting example of 
physics beyond the SM that can profoundly affect future experiments at the LHC. A hidden 
valley refers to any new sector ("v-sector") that is added to the SM with its own particles 
and interactions, neutral under the SM. A mass gap ensures that not all the particles in 
the v-sector decay to extremely-light, invisible particles. However, they can decay through 
the very weak interactions to the SM, often producing visible signals. Hidden valleys can 
lead to interesting and unexpected collider phenomenology, including hidden states decaying 
macroscopic distances away from the primary interaction vertex, or unusual final states with 
a high multiplicity of SM particles; see for example [TJ [21 E] . 

Hidden valleys have arisen in bottom-up models such as the twin-Higgs and folded-SUSY 
models [H |S] that attempt to address the hierarchy problem, and in a recent attempt to ex- 
plain the various anomalies in dark-matter searches They are also motivated by string 
compactifications: hidden sectors that are candidate hidden valleys often arise in many phe- 
nomeno logical string theory models which aim to come as close as possible to the minimal 
supersymmetric standard model, but typically have extra vector-like matter and extra gauge 
groups, see for example [TO] . 

So far various realizations of hidden valleys have been proposed, which can be generally 
categorized according to the matter content in the v-sector and the mechanism which com- 
municates the interactions between the two sectors. An interesting class of hidden valleys is 
obtained by adding to the SM a v-sector which at low energy is a pure- Yang-Mills theory 
with gauge group G v , a theory that has its own gluons ("v-gluons") and their bound states 
( "v-glueballs" ) [llj. In addition, the model contains some heavy matter fields X r , charged 
under both G v and SM gauge group, mediating the interaction between the SM and the v- 
sector. Such states are called quirks in analogy to traditional quarks^. Without additional 
massless flavors in the v-sector, the gauge group G v is confining and the dynamics of the 
v-sector will develop a mass gap, slightly larger than the confining scale A^ where the G v cou- 
pling gets strong. This scenario easily arises in many supersymmetric models with extended 
gauge groups beyond the standard model; for example, in many supersymmetric v-sectors, 

1 Quirks were considered long ago (291 130| : some of their very interesting dynamics have been studied in 

rUETJ. 
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supersymmetry breaking may lead to large masses for all matter fields. 

The spectrum of a pure Yang-Mills theory has been studied on the lattice for the case of 
a 577(3) gauge group [12]. These results are shown in figured] The spectrum of stable bound 
states includes several glueballs with different quantum numbers J PC . Some of these glueballs 
have been studied in some detail in various contexts; see for example p31 El [El EH EE3 HE1 EE] 
and a recent review [20] • The masses of the v-glueballs are determined by the confinement 
scale A„ or, equivalently, the string tension a through the relationship m ++ m 7A V m 3.7a. 
For general SU{n v ) gauge group, the leading large n v corrections to the glueball masses are 
of order l/n*. Therefore, we expect that the mass ratios in the general SU{n v ) case will not 
differ substantially from the mass ratios shown in figured] for n v — 3. All of the v-glueballs 
shown are stable against decays to the other v-glueballs, due to kinematics and/or conserved 
quantum numbers. 

The v-glueballs are non-interacting with the SM particles, except for higher dimension 
operators in the effective Lagrangian induced by the heavy mediators X r . In this work, we 
will focus on the possibility of operators between the v-sector and the Higgs sector (also 
referred to as the Higgs portal), 

0^ d (H\H)0^ (1.1) 

where H is the standard model Higgs doublet; M is a heavy mass scale, associated with the 
masses of the heavy mediator fields. Here we have split the dimension-/} operator into a 
Standard-Model part Os D of dimension D — d and a hidden-valley part Ov of dimension 
d. The (d) are constructed from gauge invariant combinations of v-gluon fields such as tr J 7 
or tr J 73 , while the oi D d ^ involve gauge invariant operators built out of the Higgs field, such 
as H'H or H^D^H . We will see that the lowest order term is given by the dimension-six 
operators of the form H^Htr J 72 . With some exceptions to be discussed below, most of the v- 
glueballs in figured] can decay via the Higgs portal interaction (jl.ip . with a strong dependence 
of the lifetimes on the confinement scale A„ and the mass scale M. 

Loops of the heavy particles also induce dimension-eight effective interactions coupling 
the v-gluons to the standard model gauge bosons, either of the form tr F 2 tr T 2 or F\ tr J 73 
where Fj (J 7 ) is the field strength tensor for standard model (hidden-valley) gauge bosons. 
The field strength tensors are contracted according to different irreducible representations of 
the Lorentz group. A detailed study of the phenomenology of these operators in the context 
of hidden valleys was carried out in [11] . In this case the v-glueball widths are dominated by 
decays into SM gauge-boson pairs, or radiative decays to another v-glueball and a photon (or 
to a lesser extent a Z boson), leading to jet and photon rich final states. 

In this paper, we will first extend the results of [llj on v-glueball decays in pure- Yang-Mills 
theory to allow for higher dimension interactions of the form (11.11) . We will need to compute 
the effective Lagrangian coupling the v-gluons to the Higgs sector. Then we will use it to find 
formulas for the decay widths of the states in the spectrum shown in figure [TJ Specifically, we 
will find that the ++ state can decay via the Higgs portal to standard model particles with 
branching ratios which are determined by the couplings of the standard model Higgs boson. 
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Other states can decay to a lighter state by emission of a Higgs boson. For v-glueball states 
which can decay either through dimension-eight operators or through the emission of a Higgs 
boson the situation is rather involved, with the relative branching ratios depending on the 
various parameters and the unknown v-glueball matrix elements. 

Our primary motivation here is in the case where the lifetimes are short enough that at 
least a few v-glueball decays can be observed at the LHC detectors. This typically requires 
the lifetimes of the hidden particles to be shorter than a few micro-seconds, if the production 
cross-section is substantial. The key point in our scenario is that with different operators 
describing the decays, often competing with each other, we typically find a large spread in the 
lifetimes of the v-glueballs, even after all the parameters have been fixed. Then there is no 
necessity to adjust any parameter to obtain short lifetimes. 

Incidentally, our results are also relevant for studies of dark matter, either in the case of self- 
interacting dark matter [21], where the v-glueballs would be the dark matter candidates, or 
indirectly, in a recent attempt to study scenarios of dark matter with novel signatures [22], 125]. 

2.5 
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1.0 
0.5 
0.0 

++ " + pc + - " 

Figure 1: Spectrum of stable glueballs in pure glue SU(3) theory [12]. Masses are shown in 
units of the lightest v-glueball mass m . 

This work is organized as follows. In section 2, we introduce our model and describe the 
effective interactions coupling the two sectors. We also classify the matrix elements. Section 
3 presents our computation of the decay modes. Section 4 is devoted to a summary of the 
different experimental constraints on the parameters of the model. Then in section 5 we present 
our numerical estimates for the branching ratios. Possible generalizations of the model are 
described in section 6. Finally, we conclude in section 7 with a brief summary of our results 
and some comments. Additional computations appear in the appendix. 
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2 The model and the effective action 



We first set up our framework and conventions. The generic scenario of hidden valleys with a 
pure-gauge hidden sector can be characterized as follows [H [H] , 

• The SM is extended by the addition of an extra SU{n v ) gauge group, with a mass scale 
scale mo in the 1 — 1000 GeV range. We will refer to this sector as the hidden valley, or 
briefly the v-sector. There are no light flavours in the v-sector, so after confinement the 
lightest states in the spectrum are bound states of v-gluons, or v-glueballs. 

• There are heavy vector-like particles X r ("mediators") that couple the v-sector very 
weakly to the visible sector at low energies. For definiteness, we will take the X r to 
transform as a fundamental representation of SU(n v ) and in complete SU(5) represen- 
tations of the standard model. We label the fields and their masses as shown in table 
[TJ Also, for convenience, we define dimensionless parameters p r = m r /M where M is an 
arbitrarily chosen mass scale, usually taken as the mass of the lightest X r particle^. 



Field 


SU(3) 


SU(2) 


17(1) 


SU(n v ) 


Mass 


*i 


3 


1 


i 

3 


n v 


m d 


Xe 


1 


2 


1 
2 


n v 


m e 




3 


1 


2 
3 


n v 


m u 




3 


2 


1 

6 


n v 


m q 


x e 


1 


1 


1 


n v 


m e 



Table 1: The new fermions X r that couple the hidden valley sector to the SM sector. 

In addition, we assume that the mediators X r can get part of their mass from electroweak 
symmetry breaking through Yukawa-like interactions. For concreteness, we consider the fol- 
lowing Lagrangian (in four- component Dirac notation), 

£mass = m rX r X r + (y^HX^ + y u X q HX u + y d X q HX d + h.c^j ^.l) 

r=d,l,u,q,e 

where H = e • . For the moment, we will restrict ourselves to the case in which the Yukawa 
couplings conserve both C and P independently. In general one or more of the Yukawa 
couplings could be CP-violating. This CP violation can contribute new terms to the effective 
action and can induce additional v-glueball decays, as will be discussed in section 6. 

There are several constraints on the couplings y r and the mass scale M from precision 
electroweak measurements. These constraints will be discussed in some detail in section 4. 

2 In this work, we normalize hypercharge as Y — Q — T3, where T3 is the third component of weak isospin. 
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We require that y r < 1 and M > 250 GeV in order to avoid potentially dangerous corrections 
to precision electroweak bounds. A high mass scale is in any case necessary to avoid current 
experimental bounds from collider searches for extra particles. 

We are concerned with the low-energy effective theory of the model described above. The 
effective interaction that couples the v-gluons and v-glueballs to the SM particles is induced 
through a loop of X particles. The coupling between the Higgs sector and the v-sector to 
non- vanishing leading order in 1/M at low energies arises from a diagram with four external 
lines, as depicted in figure [2h. This gives rise to the dimension-six operator 

2 

H6) _ a vV gtrr t r fHV (2.2) 

where a v is the SU(n v ) coupling, J 7 ^ is the v-gluon field strength tensoiEl The coefficient 
y depends on the mass ratios of the heavy particles from table [1] and their couplings to the 
Higgs field. It is given by 



y 



vt 



+ 



3^ 



+ 



(2.3) 



PlPe PqPd PqPu 

The computation of (12. 2p is presented in appendix A. Note that (I2.2p - (l2.3p are strictly valid 
in the (3 r = y 2 .v 2 H /2M 2 <C 1 limit; otherwise, the corrections to these equations can be readily 
obtained by the substitutions p r p r i — > p r p r > — (3 r in (12.31) . Once the Higgs gets an expecta- 
tion value, the following interaction terms are induced between the v-sector and the physical 
standard model Higgs boson (h), 



£(6) 



3nM 2 2 ^ 



v H htr JF lMV T^ + 



3ttM 2 

where vh = 246 GeV is the Higgs vacuum expectation valu^ 41 



(2.4) 



H 




X SM v 



SM v 




(c) 



Figure 2: Diagrams contributing to the effective action 

While dimension-six operators are only suppressed by two powers of the mass scale M, 
further suppression can arise if the operators are suppressed by small coupling y. This requires 
us to consider effective operators of dimension higher than six as well. 

3 Here we represent the v-gluon fields as = J-^JT a , where T a denote the generators of the SU(n v ) 
algebra with a common normalization tr T a T b = \8 ab . 

4 In (|2.4[) . we have not included a term which is proportional to the v-gluon kinetic term. This term can be 
removed by redefining the v-gluon coupling. 
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The next operators have D = 8, and they describe the coupling of the two v-gluons to two 
SM gauge bosons (figure [2b), either gluons (g), weak bosons (W and Z) or photons (7), as well 
as the fusion of three v-gluons into a 7 or Z (figure [2b). Within the SM, effective two gluon - 
two photon, four gluon, and three gluon - photon vertices can be found in [24] , [25] and [26] 
respectively. They were also extensively discussed in a preceding work [11]. The pure-gauge 
effective Lagrangian linking the SM sector with the v-sector reads (for our conventions see 
Appendix B) 



£(8) 



9v 



(4vr) 2 M 4 



[ g l Xl B^ B p° + 9 i X2tI pwpp° + 9 i X3tI crcr] 



X 



1 1 11 1 

^ 9pp9va + 45 e /if per + — J npQva — ~ ^p,upa 



+ 



(4tt) 2 M 4 



X 



— s^fiW - - B» v nW 

45 g P v 



• (2-5) 



Here we define Fi, 



R p2 



and F^ 



G^, which are the field tensors of the U(1)y, 
SU(2) and SU(3) SM gauge groups. We denote their couplings g iy i = 1,2, 3, while g v is the 
coupling of the new group SU{n v ). The coefficients Xi an d X encode the masses of the heavy 
particles from table [2J and their couplings to the SM gauge bosons and the Higgs doublet. 
They are summarized in table [2j The d = 4 operators (S, P, T, L) involve bilinear tr FF 
operators whereas the d = 6 operators (fl™, involves trilinears of the form tr FFF (see 
Appendix B). 





X ) Xi 


Xi 




X2 


1 , 3 


X3 


1,1,2 
^ + ^ + ^ 


X 


1 1 2,1,1 
^-^f-^ + ^ + ^J 



Table 2: The coefficients % arise from a sum over the SM charges of X particles running in the loop. 
The Xi, i = 1,2, 3, arise from the diagram in figure^b) with two external SM gauge bosons of group 
i, while x 1S determined by the diagram [2](c) with a single hypercharge-boson on an external line. 
Here p r = m r /M. 

We will not consider operators of dimension 8 containing both gauge-boson fields and the 
Higgs field, as well as dimension-ten or higher operators, since their effects are suppressed by 
extra powers of M and/or y. 

The interactions in the effective action then allow the v-glueballs that cannot decay within 
the v-sector to decay to final states with SM particles and at most one v-glueball. This is 
similar to the weak decays of hadrons, such as n + — > Iv, n — > pu e e~ , and so forth, where the 
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Fermi effective interaction allows otherwise stable hadrons to decay into leptons. To compute 
these decays, we will only need the following factorized matrix element^]: 

(SM\Oi D -V\0)(0\Olf)\e K ) , (2.6) 

(^M|O( D - d )|O)(0 K ,|O^|e K ). (2.7) 

Here d is the mass dimension of the operator in the v-sector, (SM\ schematically represents 
a state built from Standard Model particles, and \Q K ) and \Q K ') refer to v-glueball states 
with quantum numbers k, which include spin J, parity P and charge-conjugation C . The 
SM part {SM\O { s D ' d) \0) can be evaluated by the usual perturbative methods of quantum field 
theory, but a computation of the hidden-sector matrix elements (0\Oi |0«) and (Q K '\Oi, \® K ) 
requires the use of non-perturbative methods. 



2.1 Matrix elements 

We wish to classify the non-vanishing v-sector matrix elements of the scalar operator S = 
tr J z llv J r ^ v in (12.2ft . As we saw, the matrix elements relevant to v-glueball transitions are given 
by (O|5|0 K ) and (Q K r\S\Q K ), where K and Q K > refer to v-glueball states with given quantum 
numbers. 

It is convenient to write the most general possible matrix elements in terms of a few Lorentz 
invariant amplitudes or form factors. The decomposition of S into irreducible representations 
of the Lorentz group contains only the ++ quantum numbers [28] . This allows the ++ state 
to decay directly to standard model particles. For the annihilation matrix element of the ++ 
v-glueball, we will write 

(0|5|0 ++ )=F^ + (2.8) 

where is the ++ decay constant. 

Likewise, we can consider the transition matrix elements of the operator S between a spin 
J state with momentum p and a spin J' state with momentum q. For the moment we make 
no assumptions about the parity or charge conjugation quantum numbers of the states. Of 
course, if we impose parity some of the transitions may be forbidden. The matrix elements 
can be compactly written as 

(J'\S\J) ee J2 M fj' M *J> (2.9) 

i 

where now Mjj' is the transition matrix, which depends on the transferred momentum, and 
M-jji is determined by the Lorentz representations of | J) and | J'). In Appendix C we have 
listed Ai jji for the simplest cases considered later in this work. 

5 As mentioned in the introduction, decays with no SM particles in the final state are forbidden because of 
kinematics and/or conserved quantum numbers. 



The main uncertainties in the study of decays of the v-glueballs stem from the evaluation 
of the transition matrix elements Mjj 1 ? which at present are unknown. Many of these could be 
in principle be determined by additional lattice computations. In spite of these uncertainties, 
we will still be able to obtain many interesting and robust results. 

3 Decay rates 

An interesting feature of the model described above is that it has a spectrum of v-glueballs 
which can interact very weakly with the particles in the standard model through effective 
operators of different mass dimensions. The dimension-six operator, as we see below, permits 
the v-glueballs to decay directly into standard model particles (the ++ ) or radiatively by 
emitting a Higgs boson (the 2++, 2~+ 3++, 3+~, 2+~, 0+" 1"", 2", 3—). On the other 
hand, dimension-eight operators contribute with additional decay modes. As shown in [TT] . 
these include direct annihilations into SM gauge-boson pairs (the ++ , 2 ++ , h and 2~ + ) or 
C-changing radiative transitions with emission of a photon or a Z boson (all others). 

In this section, we will compute the decay rates for most of the v-glueballs in figure [TJ 
We will first present our results for the v-glueball decays induced by dimension-six operators 
given in (12. 2p . Then we will make a quick review of the results of [TT] for dimension-eight 
operators that are relevant for our work. 

3.1 V-glueball decays by dimension-six operators 

We begin with the ++ v-glueball, which can be created by the operator S. The dimension-six 
coupling leads to ++ — > (( annihilations, where ( denote any of the final states of the Higgs 
boson. Above the threshold for Higgs boson pair production, the ++ — > hh decay can also 
proceed with a sizeable rate. Other decay modes induced by dimension-six couplings include 
processes of the form 6 K — > B K >h, where K , K / denote two v-glueballs with given quantum 
numbers. As an example of a Higgs-radiative decay, we will present the computation of the 
2 ++ — > ++ h decay with some detail. Then we will consider the general case K — > Q K >h for 
v-glueballs with arbitrary quantum numbers. In the end, we make some comments on the h 
and 1 H v-glueballs, which are the only ones that are not permitted to decay via dimension-six 
operators. 

Annihilation of the ++ state 

The scalar state can be created or destroyed by the S operator. Then, the effective interaction 
(12. 2 p allows the decay of the ++ state via s-channel Higgs-boson exchange ++ — > h* — > ((, 
where ( collectively denotes a standard model particle. According to (12. 2p . the amplitude for 
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this decay reads 



(CCK/7 + m%ZpZ» + 2m 2 w W^WnO) ^r^| (Ol^lO (3-1) 
where a v = g 2 /(4:7r) and m# is the Higgs mass. The width of the decay is given by 

(2 \ 2 

where F^ + = (0|S'|0 ++ ) is the ++ decay constant. Here TfJ^^{m 2 + ) is the width for the decay 
h — > (( for a standard model Higgs boson with a mass m + . Then, formula (I3.2p implies that 
the branching ratios of the ++ are those of the SM Higgs boson in the range of mass of 
interest. Expressions for the branching ratios and full width of the Higgs boson can be found 
in the literature, for masses ranging from a few MeV up to 1 TeV (for a review, see |27j). 

Although these are standard model tree-level results, we should also remark that the 
range of validity of f)3.2p is beyond the simple perturbative QCD domain. For masses below 
2 — 3 GeV, the ++ v-glueball can decay into a pair of hadrons via its interaction with two 
gluons through a top-quark loop or its interaction with quarks. The hadronization of these 
quarks and gluons is a rather complex and non-perturbative process. However, the branching 
ratios are still given by (13. 2p . with ( now running over the possible hadrons in the final states, 
such as 7r, K, and so forth. 

If m + > 2m^, the decay channel ++ — > hh opens up and the partial width is given by 

1 fy 2 a v F^ + \ 2 f 3m 2 H \\, 2 2 2 yii/2 to o\ 

where g(x,y; z) = (1 — x/z — y/z) 2 —Axy/z 2 . The off-shell decays ++ — > h*h and + — > h*h* 
are also possible in the intermediate mass range. However, these channels receive an extra 
suppression from the smaller available phase space, so it is reasonable to expect they will have 
little effect on the full ++ width. 

As m + becomes larger than 2m t and 2m#, we have the following approximate relationship 
among the dominant decay rates 



T w+W - : Y Z z :T hh :T t i = 2:l:l:?>Xt \?l - x t (3.4) 

where x t = ^m 2 / m 2 + . 



An example of Higgs-radiative decays: the 2 ++ — > ++ h case 

As an illustration of the computation of the decay rates for transitions with the emission of 
Higgs boson, let us first consider in some detail the decay of the 2 ++ v-glueball. Contrarily 
to the ++ , direct annihilation of the 2 ++ v-glueball into Higgs bosons would require an 
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operator in the effective action of dimension D = 8, and is hence negligible. Instead, for 
m 2 + — m + > m H , the S operator induces the 2 ++ — > ++ h decay. In a parton-model picture, 
one may imagine this process as being caused by the radiative emission of a Higgs boson 
through g v — >■ g v h, with one spectator v-gluon going to the final state. The amplitude of this 
two-body decay is given by 

|^< % |0><0~|S|2~) ^* !(V(Y (3.5) 

where g M is the momentum of the ++ and the polarization of 2 ++ . Here Mq +2+ denotes 
the transition matrix element defined by Mg +2+ e Ml/ g /i g l/ = (0 ++ |S'|2 ++ ); for simplicity we have 
assumed it is independent of the transferred momentum. The corresponding width of the 
decay reads 

1 = 48o/(3^M^ + [ 5 K + ,<;^)] 5/2 (M s +2+ ) 2 (3.6) 

If the 2 ++ is not heavy enough, m 2 + —m + < rrifj, it can decay to the ++ and SM particles 
through the emission of an off-shell Higgs. The corresponding amplitude for the three-body 
decay 2++ ->- 0++CC is 

^ (CCK/7 + niZ^ + 2<W-W- 1 0) W —^L_ ¥T _ (0 + + | S |2+ +) (3.7) 
where k 2 is the transferred momentum. The width is given by 



,3 



M s +2+ ) 2 / V a v v H \ 2 r 2 [j^.^iO] 5 ' 2 ^ ( ffll2 ) 



_ m 2+l ivi 0+2+ ; / y ut v u H \ I , 2 \P V '"0+ ' " "12 ) " "2+ ) J ^ h^CC^"" 1 ^ ("? «1 
1 2++ ^0 ++ CC - ^ I^^M^ J 7 ^ ( m 2 2 _ m ^ )2 + (r 5M )2m ^ ^ ' ' 

where the limits of integration are set by available phase space. 
Higgs-radiative decays. General case J — > J'h 

We now consider the 2-body decays 0j — > Qj/ h, in which Bj, Qjr are v-glueballs with spin 
J, J' respectively. For the moment, we make no reference to the parity of the v-glueballs 
and proceed generally. Parity conservation will force some of the matrix elements to be zero. 
Making use of the general formulas for the matrix elements (1A-23l) - (1A-30j) . we arrive at the 
decay rate 



r?U«) = ^Ji'tlor^ r^ M U 2 r% [g{m^m%;m^ (3.9) 



16vr(37r) 2 M 4 mj (2 J + 1) 



where i runs over the various form factors. The coefficients Tjj, are dimensionless functions 
of the masses and depend on the angular momentum transfer associated with each transition. 
They are summarized in table [31 
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1 

1 


2 


3 


4 


5 


1 32 


^(4x z + 28x + 35) 


f^(x + 2) 


~L5~ 


f^(2x + 7) 


4x 4 
15 


r 31 


2x^ / O i T\ 

f^(3x + 7) 


2x 4 
5 


8x 3 
~L5~ 






p(») 
1 30 


2x 3 
5 










r (i) 
1 22 


1 / a 9 i on i An\ 

^(4x + oUx + 45) 


Yg(8x + 17) 


4x 4 
9 


f (x + 5) 


X 3 


p(i) 
r 21 


3 (2x + 5) 


2x 3 
3 


™2 
X 






r (i) 
1 20 


2x 2 
3 










p(i) 
1 11 


x + 3 




2x 






r w 

1 10 













Table 3: The coefficients T jj, arise from the average of the squared matrix elements. We denote 
x = ■^-g(m 2 j,,m 2 H ;mj). The dashes denote those cases where form factors are absent. 

Below the threshold for Higgs boson production, the decay rate for the 3-body decay 
Oj Qji CC reads, 

rj-fj'« = - / dml 2 m 12 r%j, h (m 12 ) rff cc (m 12 ). (3.10) 

71 J / -H'''l2> " L H) 

where A(m\ 2l m 2 H ) = {m\ 2 — m 2 H ) 2 + m 2 H (T^ M ) 2 . The integration automatically includes the 
case where the radiated Higgs boson can be close to onshell. 

Some of the transitions may be parity-forbidden (e.g. 0~ + — > V h). However, it is a 
straightforward exercise to check that for each v-glueball G K in figure [1] , there exists at least 
one other v-glueball Q K / such that the transition K — > Q K >h is allowed, with a rate which is 
given by f l3.9p - fl3.10p . The only exceptions are the +_ and l + ~ v-glueballs that we discuss 
next. 

Decays of the pseudoscalar and pseudovector 

The only states that are not allowed to decay via dimension-six operators are the l + ~ and the 
0~ + v-glueballs. 

Since the l + ~ state is the lightest state in the C-odd sector, it necessarily has to decay to 
a v-glueball of opposite C. One possibility is that the l + ~ decays by radiatively emitting a 
Higgs boson (e.g. l + ~ — > 0~ + h). However, this decay mode would violate C and, hence, is 
forbidden. 

A second way to induce decays of the 1 H v-glueball would be via its coupling to the 
hypercharge current H^D^H. In this case, there are three v-sector operators that can be 
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contracted with H^D^H, namely, 

trF^D^, \xJ^D p T alt , \zTtfD a F#. (3.11) 

However, one can see that these operators cannot induce C-changing transitions. First, notice 
that classically tr T^D^T 01 ^ = fc^tr T^T 01 ^ . As explained in [23], this implies that the 
transitions induced by the first operator in (13. lip are not new. They are just the same ones 
created by the operator S. Likewise, using equations of motion and conservation of the energy- 
momentum tensor, the second operator in (13. lip can be related to a total derivative of the 
operator S. Finally, equations of motion also imply that the last operator in (13. lip vanishes 
identically. Therefore, up to operators of total mass dimension six, the l +_ state is stable. In 
the next subsection, we shall see that dimension-eight operators can induce photon-radiative 
decays of the l + ~ v-glueball to C-even v-glueballs. 

On the other hand, C-invariance by itself would allow the Higgs-radiative transition 0~ + — > 
++ h. However, this decay could not conserve both angular momentum and parity P: since 
the initial state is 0~, angular momentum conservation requires the orbital angular momentum 
of the + and h final state to be L = 0, which in turn requires total parity P = +1, rather 
than P = — 1 as demanded by parity conservation. This decay mode is thus forbidden. This 
is analogous to the way that the rj meson strong interaction mode rj — > tttt is forbidden in the 
SM. A similar argument shows that the three-body decay h — > ++ hh is not permitted. In 
this case, the corresponding rj — > tttttt decay is allowed in the SM because, contrary to the ++ 
and h, the pions have intrinsic parity —1. As we will see in section 6, this line of argument 
alters if we relax our assumptions and allow for P-violating couplings. 

3.2 V-glueball decays by dimension-eight operators 

The operators in the effective action (12.51) induce decays for all of the v-glueball states in figure 
[TJ The lightest states in the C-even sector (the ++ , 2 ++ , 0~ + and 2 h ) can directly decay 
to pairs of standard model gauge bosons (gg, 77, ZZ, WW or 7Z) via the S, P, T and L 
operators. All other states can decay by radiatively emitting a photon or, to a lesser extent, 
a Z boson, via the d = 6 D = 8 operators fl^ v . Here we briefly summarize the results of [UJ 
concerning the computation of the decay rates induced by D = 8 operators. 

In order to retain simplicity, we will often assume, in the subsequent discussion, the X 
fields form approximately degenerate multiplets of SU(5), i.e. p r ~ 1. When there is a large 
hierarchy between the colored and the uncolored X particles, Pd^> Pi = Pe and p a ~ p q ^> p e , 
the decay pattern becomes slightly more complicated because the decay channels into 77, 7Z, 
ZZ and WW may all play a role and even dominate in some regions of the parameter space. 
We will make some comments later in this paper on some of the interesting phenomenology 
that arises in this regime. 

For the decay of the ++ , 2 ++ , 0~ + and 2~ + v-glueballs into gluons we have the following 
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rates, 



r < 0+ -^ = l§4s) 2m ° +(F ° ++)2 ' < 312 > 



r < 2+ - 99 > = 



i niv ,_ T , 2 4/1 



2 \4SI > F2++) + 3 IsO 1 (F2+ " 



2 



(3.13) 



t ^^^ = M^(ie) < (f °-+> 2 (3 ' 14) 
r ( 2 --^^^-(s) 2(F -) 2 (3 ' 15) 

where a s = g^/(4ir) is the QCD coupling constant and the coefficient X3 is given in table [2j 
Of great interest is the branching fraction into two photons, 

r(e^77) = ia 2 xl, f316 , 

T(Q^gg) 2a* sX l 

for 6 = ++ , 2 ++ , CT + , 2~ + . Here a is the fine structure constant and x-y = Xi + X2/2, where 
Xi and X2 are shown in table EJ The expressions for the branching fractions into electroweak 
bosons are omitted for the sake of brevity but can be found in [TT]. Some comments on the 
weak boson decay modes will be given in section 5. 

In the C-odd sector, all the v-glueballs decay radiatively with the emission of a photon to 
the lightest v-glueballs in the C-even sector. Direct annihilations into three SM gauge bosons 
are suppressed since they would be induced by dimension-twelve operators. One can also show 
that three-body decays are suppressed due to the small phase space that is available for these 
decays [TT] . 

The lightest states in the C-odd sector are the pseudovector l + ~ and the vector 1 . The 
width of the decay l + ~ — > ++ + 7 is given by 



aa l 2 ( m i+ _ m2 



i+_>o++ 7 = n , " o X 1 3 ° + ( M "-o++)^- (3-17) 



24ttM 8 



m 



A similar expression holds for the decay 1 — > ++ + 7 



i,3 ( iy-n 2 



^ v 2 V" 1 !- nb Q+) f-xjSi \2 /n lfi \ 

ri ^ 0++ ^ " 24^ X ^ ( M i-o ++ ) • (3-18) 



For the 1 state, annihilation to SM fermions via an off-shell 7 or Z is also possible, with 
branching ratio 

r = (Yl + n) 3 (_^_) 2 . ,3,9) 

tl-^0++7 COS Vw \ m i- — 771 0+ / \ iV1 l--0++/ 
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Here Yl and Yr are left and right hypercharge of the emitted fermions. Decay to electrons 
and muons will be reconstructable as a resonance, so despite its uncertain branching fractions, 
this decay mode is mode is worthy of careful consideration. 

One can also generalize formulas ( 13.171) and ( 13.181) to include the radiative decays of the 
heavier C-odd v-glueballs. Easy computations show that 

1 J ^ ++ " ~ 4vrM8 2^J(2J)!(2J+l) m 2 / +1 m^r 2 { J0++) ( } 

aal 2 K-4T 2 J -\J\f 
J ^ 2++1 48ttM8 X m 2 / +3 ml J + 3 J (2 J)! (2 J + 1) 

(2(71 J + 65)m^ ++ + 3(5J + 3)m^ + 3(5J + 3)m^ ++ ) (IVr? 2++ ) 2 . (3.21) 

with similar expressions for the modes J — > 0~ + 7 and J — > 2~ +/ -f. 



3.3 Summary of decays 

In table H] we summarize the final states for the most important decay channels of the v- 
glueballs in figure [1] for D = 6 and D = 8 operators. We therefore see the presence of 



State 


D = 6 operators 


D = 8 operators 


0++ 


bb, W+W-, ZZ, hh 


gg, WW, ZZ, Z 7 , 77 


2 ±+ 


±+ h(h*) 


gg, WW, ZZ, Z 1 , 77 


o-+ 




gg, WW, ZZ, Z 1 , 77 


3++ 


0~+h, 2 ±+ h{h*) 


0- + gg 2 ++ gg, l + ~7 


1+- 




±+ 7, 2-+7 


1— 


l+-h(h*) 


±+ 7 , 2 ±+ 7 , // 


0+-, 2+-, 3 + - 
2-, 3- 


J p -h(h*) 


±+ 7, 2 ±+ 7 



Table 4: Possible final states of the various v-glueballs in figure [T] generated by D = 6 and 
D = 8 operators. Note the absence of Higgs- mediated decay modes for the _+ and l + ~ 
v-glueballs. Here J p ~ denotes a C-odd v-glueball state. 

operators of different mass dimensions opens a plethora of decay modes, which is particularly 
interesting from the phenomenological point of view, but complex to analyze. It is the purpose 
of the next two sections to disentangle the effects from D = 6 and D = 8 operators, and extract 
the most frequent decay modes. 
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4 Constraints on new physics 



In this section, we discuss direct experimental constraints on the operators used in this paper 
from Tevatron searches for new physics, as well as potential limits from precision electroweak 
measurements. 

4.1 Direct searches 

In order to translate nonobservation of X particle events above the standard model back- 
ground into bounds on the X mass, it is important to understand the possible collider signals 
associated with the production of X particles at hadron colliders such as the Tevatron or the 
LHC. As was pointed out in [291 EQ], and more recently in [TJ I3T] . theories with quirks like 
the X particles give rise to a quite unusual dynamics. Here we briefly rephrase and elaborate 
some of the results of [TJ [31] to agree with the regime we are considering in this paper. 

At the Tevatron, color-singlet X particles could be pair produced through an off-shell 
photon, Z, or W. If the X particles are colored, their pair production through an off-shell 
gluon is also possible. To a great extent, the production rate is only dependent on the gauge 
couplings of the X particles and so is fixed by their standard model quantum numbers. The 
X and X are initially produced with considerable kinetic energy. Due to v-confmement, they 
form a bound state of X-onium with an SU{n v ) string (flux tube) being stretched between 
the X and the X. As the X and X move apart, the potential energy stored in the string 
increases, but contrary to QCD strings, the SU(n v ) string cannot break, since there are no 
light v-quarks in the v-sector whereby the string can split. Instead, the X and X oscillate back 
and forth losing their kinetic energy and angular momentum by radiating some combination 
of photons, hadrons and v-glueballs. Eventually the X and X recombine and annihilate back 
into lighter states. 

In the case of colored X particles, the annihilation is most often into v-gluons, which at 
long-distance become two or more v-glueballs. As explained in section 3, the v-glueballs can 
then decay via a loop of X particles, producing photons, gluons, bottom quarks, and so forth. 
As will be seen in section 5, v-glueball decays typically result in an observable signal, providing 
A-u ^ 1 GeV. For color-singlet X particles that are produced via an off-shell photon or Z, 
the primary annihilation channel is also to v-glueballs, which can decay back to SM particles. 
However, the dominant cross-section for production of color-singlet X particles is through an 
off-shell W boson [8]. The resulting X + X° bound state is electrically charged and cannot 
annihilate into v-glueballs alon^j. Instead, it may annihilate to leptons or quarks via an off- 
shell W, to Wj or WZ, or to Wg v g v , that is W plus some v-glueballs. Other contributions to 
the final states arise from the v-glueballs radiated during the X-onium relaxation [TJ |8j [3Tj . 
These might help increase the discovery reach and deserve further investigation which we leave 
for future work. 

6 One can show that the splitting between X + and X° is so small that X + is stable against decay to X° 
on the time scales of interest; for example, see |S]. 
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The W"f process might allow X-onium to be discovered directly. Events with a photon, 
lepton and missing energy ( I + 7 + fir) would reveal a X-onium resonance in the transverse 
mass distribution above the W~f background. On the other hand, when the v-glueballs decay, 
they may produce spectacular collider signals, such as pp — >■ XX — >■ 7777, 77r + r~. The 47 
channel, in particular, is very clean, since it is essentially background free, and would reveal 
the v-glueballs, if the rate is large enough to be detected. 

There have not been specific searches for quirks as yet. Present limits are in general very 
mild because the cross-section for production of quirks at the Tevatron becomes very small 
once the quirk masses are larger than ~ 250 GeV and because not all the collected data has 
been analyzed yet. Below we present bounds on the masses of the X particles from CDF 
searches for new physics in I + 7 + JS? events and in searches for anomalous production of two 
photons and at least one more energetic, isolated and well-identified object (r lepton or 7). 
This analysis by no means is exhaustive, but is meant to illustrate potential discovery signals 
and the bounds that can arise. 

Our estimates below are based on the following assumptions. First, we determine the 
rate for production of X particles through perturbative computations; the cross-section a xx 
is estimated as a function of the mass M x and for n v = 2, both for the cases of colored 
and uncolored particles. Second, motivated by a thermal model of fragmentation [22], we 
have estimated the relative probability of producing any v-glueball in X-onium annihilations 
assuming production of the ++ at 50 — 60%, and that of the 2 ++ , _+ and 2 h totaling 
40 — 50% in all. Finally, for the purpose of obtaining an approximate bound we assume a 
100% efficiency for detecting v-glueballs. The last two assumptions are unlikely to be accurate, 
as we will discuss more thoroughly in our LHC study, but they are intended to give a rough 
estimate for the bounds. 

(i) l + 'j + JSx'- The CDF collaboration has searched for the anomalous production of events 
containing a high-transverse momentum charged lepton and photon, accompanied by missing 
transverse energy [32]. Using 929pb _1 of CDF Run II data, these searches exclude X masses 
below 200 GeV for color-singlet X particles and 250 GeV for colored X particles, under the 
assumption that the X-onium has lost most of its energy by radiation of photons and/or 
v-glueballs and the annihilation takes place at or near the ground state. 

(ii) 27+7 and 27 + r: The CDF collaboration has also searched for the inclusive production 
of diphoton events [33]. Using 1155 pb -1 of integrated luminosity, the diphoton plus third 
photon search places lower limits on the X mass below 200 GeV for color-singlet X particles 
and y ~ 0. The corresponding limit on the X mass for colored X particles turns out to be 
less severe, roughly below 150 GeV, because, in this case, the branching fraction to photons 
is < 0.35%. Similarly, using 2014 pb _1 , the diphoton plus tau search places lower limits on 
the X mass for y ~ 1 below 250 GeV for color-singlet X particles and 175 GeV for colored 
particles. 

Some comments are in order. The mass bounds from potential v-glueball signals are 
necessarily model dependent. Unfortunately, unlike our computations of the branching ratios, 
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there is no way to make a reliable estimate of the probability for producing any given v-glueball 
state in X-onium annihilation, since phenomenological input from QCD is not relevant to the 
pure-glue case. V-glueballs could be created during the relaxation of X-oniums, and during 
their annihilation, so events with more than two v-glueballs might be common. Generally, 
with more than two v-glueballs produced in each X-onium annihilation, events with multiple 
photons would be common and easier to observe. However, the events may be more cluttered 
and the energy of the photons lower, leading potentially to lower detection efficiencies. Since 
we cannot model this reliably, in the simple analysis of this paper, we will rely on the model 
independent bounds from the production of X-onium through s-channel outlined above. 

Moreover, the bounds from the W+ photon final state assume that the X-onium state 
decays to its ground state before annihilating. However, it has been argued in [31 J that the 
annihilation may often occur at higher energy, before the ground state is reached. With 
many X-onium states annihilating at different energies, the signal in the transverse mass 
distribution will be diluted. Besides, annihilation to two fermions via an off-shell W or to W+ 
v-glueballs becomes increasingly significant at higher energies. While the branching ratios to 
Wj/Z, f f and W+ v-glueballs can be estimated, the wide range of X-onium states and the 
uncertain annihilation probability as a function of the energy make any precise evaluation 
almost impossible. Therefore, the Tevatron limits are in fact weaker than suggested above. 

In passing, we note that the bounds on a fourth generation of particles from Higgs searches 
at the Tevatron do not apply to heavy fermions that get only part of their mass from elec- 
troweak symmetry breaking, such as the X particles in our paper. For an additional pair of 
fermions that get most of their mass from electroweak symmetry breaking, the Higgs produc- 
tion cross-section & gg ^H is known to increase by a factor of roughly 9, making a dramatic 
effect on Higgs physics [33]. This allowed CDF and DO to rule out fourth generation quarks 
for a Higgs mass in the window 145 GeV < m# < 185 GeV. However, for heavy particles 
that get only part of their mass by electroweak symmetry breaking, their contribution to the 
cross-section <J gg ^,H decreases and the bounds are generally much weaker or absent. 

Before moving on to the discussion of electroweak precision constraints we should briefly 
comment on the possibility of producing v-glueballs in Higgs decays. Up to now we have 
considered the production of v-glueballs as a by-product of X-onium annihilation and relax- 
ation. However, there is another possibility for this production to occur. For sufficiently small 
v-glueball masses, the interaction (I2.2p can also mediate processes such as h — > K K and h — > 
O K O K K . Assuming only minor phase space suppression, the branching ratio to v-glueballs for 
a SM Higgs below 140 GeV and n v = 2 - 4 is of order (a v y 2 vlm h /37rM 2 m b ) 2 ~ 10~ 3 - 10~ 2 . 
This implies a cross-section for producing v-glueballs at the Tevatron of order a few fb. Even 
though the cross-section is small, when combined with X-onium annihilations, the production 
of v-glueballs in Higgs decays may be useful in the low mass range m — 1 — 70 GeV. The pro- 
duced v-glueballs can then decay to a wide variety of final states, including b quarks, r leptons, 
and multiple photons, among others; for small v-glueball masses, a significant fraction of these 
decays may occur with displaced vertices. At the Tevatron, there may be potential to observe 
a few events in the bbr + r~ or bbbb channels, although with limited statistics. In the case of 
displaced vertices, the DO collaboration has searched for pair production of neutral long-lived 
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particles decaying to a bb pair, and found no significant excess above the background [55] . 
This search could only exclude branching fractions of order 1, so that in our case it does not 
imply any bounds. At the LHC, the background problem is more severe, but with a larger 
cross-section and a larger integrated luminosity than at the Tevatron, the observation of a few 
v-glueball events may also be possible. Depending on kinematics, the case of photons that 
originate away from the primary interaction point may also play an interesting role. A search 
for displaced photons has been done at DO [37], and may be possible at ATLAS, due to the 
longitudinal structure of the electromagnetic calorimeter (ECAL). Besides, CDF and CMS 
allow the identification of delayed photons from long-lived particles using ECAL timing infor- 
mation (see [31] for results obtained using the CDF detector). The question then is whether 
at a rate ~ 1CT 3 — 1CT 2 the LHC will have some sensitivity. Keeping these non-standard Higgs 
decays in mind we leave the analysis of the discovery reach for future work. 



4.2 Electroweak oblique corrections 

If new heavy particles exist, they can manifest in the standard model in terms of corrections 
to the gauge-boson self-energies. When the new physics scale is much larger than Mz, this 
effect can be described by just three parameters (S, T, U) at the one-loop level [3SJ EH] [J 

S = 16vr^ [n 33 (g 2 ) - n 3Q (g 2 )] 1,2=0, (4.1) 

47T 

T =7^[ n ii( )- n 33(0)], (4.2) 

U = 167T A [n n (g 2 ) - n 33 (g 2 )] | ?2=0 , (4.3) 

where My/ is the mass of the W, and s 2 = sin 2 6w- The subscripts 1 and 3 refer to the weak 
SU(2) currents, while Q denotes the electromagnetic current. In practise, only the S and 
T parameters are relevant for our work because the U parameter is suppressed by an extra 
factor of the heavy fermion masses. These parameters are a measure of the size of electroweak 
breaking, which is parametrized by the breaking scale yvjj- On the other hand, they must be 
suppressed by the mass scale M in the limit M — » oo. Therefore even if one introduces heavy 
fermions with extremely large couplings, their contributions to the S, T and U parameters 
can become small at least in the limit M 3> yvu- This leaves ample available parameter space 
within which extra vector-like fermions are in agreement with all experimental constraints. 

In this section, we analyze how current limits on S and T from precision electroweak fits 
can be used to obtain constraints on the mass splittings of heavy fermions. To calculate S 
and T, we use exact one-loop expressions for the gauge-bosons self-energies which are valid 
for all values of the new vector-like fermion masses. 

7 The reader should not confuse these S and T which are both vacuum polarization functions with the 
operators S and T in (|2.5|) 
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Our scenario contains three vector-like fermion fields, one doublet and two singlets of 
SU(2) L , transforming under the fundamental representation of SU(n v ), 

^ = ( % ) ^u = U ipd = D (4.4) 



. Q 

with hypercharge Y, Y + 1/2 and Y — 1/2, respectively. We will make our computations using 
the above quark-like fermions. With a simple modification, our results can be readily applied 
to the case in which the fermions have lepton quantum numbers. 

The full fermionic Lagrangian is given by 

£ = i) q {D^ + m q )ip q + ipuiD^Y + m u )^ u + MD^Y + m d )ip d + 

+ (yui> q Htp u + y d i> q H^ d + h.c.) (4.5) 

where m q , m u and m d are Dirac masses. For simplicity we have assumed that the mass matrix 
is symmetric and real, but complex masses may be present as well. Finally, the covariant 
derivative is 

= 8, - igT a W; - ig'YB,. (4.6) 

When the Higgs acquires an expectation value (H) = ^ ^ ^ , off-diagonal mass terms are 
induced for ip q , ip u and ip d 

M u =( m « VuV ) M d =( m « VdV ) (4.7) 
V VuV m u J \VdV m d J 

The mass matrix M is diagonalized by 

( tl ) = ( -Sl c! ) ( X ) ( tl ) = ( -S2 C2 ) ( X ) (4 ' 8) 

where C\ = cos0i, Si = sin0i, c 2 = cos0 2 , s 2 = sin0 2 with 

^VuV , 2y d v 

tan 20i = — tan 202 = — • 4.9 

m q -m u m q -m d 

The corresponding eigenvalues are given by 

"ii,2 = ^ (rn q + m u ± yj (m, - m u ) 2 + 4y%J (4.10) 

™3,4 = ^ (m q + m d ± (m q - m d ) 2 + (4.11) 
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with the following inverse relations, 



m q = c-l mi + s-l m2 = c 2 m 3 + s 2 m 4 (4-12) 

m u = s\m l + c\ m 2 

m d = s 2 m 3 + c 2 m 4 

2y u v = (mi - m 2 ) sin20i 

2y d v = (m 3 - m 4 ) sin20 2 . 

These expressions now permit the evaluation of the S and T parameters as follows, 

? = ^ {- C j(Y + |) log-? - ^(Y + f ) logm^ + |(y - |) logm^ 

+|(Y - |) logm 2 + Mfi n'(0, mx, ma) + ^n'(0, m 3 , m 4 ) | (4.13) 



n iV 

T = 8ttJ 2 M 2 ( c i c 2 n ( ' m i' m 3) + CiS 2 II(0, mi, m 4 ) + SiS 2 IT(0, m 2 , m 4 ) + 

+SiC 2 Il(0, m 2 , m 3 ) — CiSiIl(0, mi,m 2 ) — c 2 s 2 II(0, m 3 , m 4 )} (4.14) 
where the functions 11(0, mi, m 2 ) and n'(0,mi,m 2 ) are given by 

II (0, mi, m 2 ) = — ^ k (— mf + 4m?m 2 + 4m? (mi — 2m 2 ) logmi — 4mim 2 + 

mi — m?> v 



+4m 2 (2mi - m 2 ) logm 2 + m 2 ) (4.15) 



II' (0, mi, m 2 ) = — — t. 7—r (-12m?(m? - 3mim 2 + 3m 2 ) logmi+ 

9(mf — m 2 ) d v 

+12m 2 (3m^ — 3m?m 2 + m 2 ) logm 2 + 

+ {m\ - ml)(2m\ + 9mfm 2 - 16m 2 m 2 + 9mim 2 + 2m 2 )) . (4.16) 

In this context, the S and T parameters are a measure of the deviation of the heavy 
particles from the pure Dirac mass case. When the couplings y r are turned off, the custodial 
SU(2) C and isospin SU(2) L symmetries are restored, and both S and T vanish. To see this 
in more detail, we can expand the S and T parameters in powers of (m 1 — m 2 )/(mi + m 2 ), 
(m 3 — m 4 )/ (m 3 + m 4 ). If m q = m u = m d = M ^> y u v, y d v we obtain, 

N c n v v 2 [(11 + 20Y)yl + (11 - 20Y)yj] 
30vrM 2 

N c n v v 4 (yl-y 2 d ) 2 , ■ 
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Figure 3: The bounds at 95% CL on the (M,y) parameters from constraints on the oblique 
parameters (S,T) for n v = 2,3,4 and two different regimes: p r ~ 1 (y = ylyi) (Left panel), 
and pj = p n = p q 3> Pi = p e ~ l(y = Hi) (Right panel). The upper-left region is excluded in 
these plots. 

where sw = sin 9w and My/ is the mass of the W boson. Then we see that S, T — >■ when 
y u ,yd — > 0, i.e., m 2 — > mi and m 4 — > m 3 . The corresponding formulas for the lepton-like 
fermion case can be obtained by substituting y u — > 0, yd — > yi, N c — > 1 in equations (I4.17P 
and dUHD - 

Fits of the combined electroweak data provide constraints on the S and T parameters and 
have been obtained in many places. Here we use the results from the PDG fits. The standard 
model is defined by (S, T) = (0, 0) with m t = 170.9 GeV and mu = 115 GeV. The best fit to 
data is (without fixing U = 0) 

S = -0.10 ±0.10 (4.19) 

T = -0.08 ±0.11. (4.20) 

These equations then imply the bounds S < 0.10 and T < 0.13 at 95% CL. From the ex- 
perimental data, we can readily find the constraints on the parameters of the model. For 
example, for lepton-like fermions, if we take Y = —1/2 and mi = 1 TeV, then S ~ O.Ol^ 2 ?^ 
and T ~ 0.08yfn v . Notice that the small value of S leads to only mild constraints on yi and 
n v . By constrast, the parameter T provides the most stringent bounds. 

There is, however, one important effect that can mitigate the contribution to T. From 
(14. 181) , we see that electroweak corrections to the T parameter are insensitive along the y u = yd 
direction in the parameter space, where the isospin symmetry is restored. In this case, the 
only constraints on the y u and yd contributions to y are given by the S parameter, which, as 
we have just seen, are less stringent. 

In figure [3] we show the bounds on y at the 95% CL. from the X particles in table [1] 
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as a function of the mass scale M and for n v = 2,3,4. For illustration, we have considered 
Vu = Vd — Vi an d two different regimes: the degenerate case, p r ~ 1, (Left panel) and the 
non-degenerate case, pj = p„ = p q ^> pi = p e w 1, (Right panel). In this case, one has 
1/ = a/7?// and y = yi for the degenerate and non-degenerate cases, respectively, (cf. (12. 3p ). 
Since y u and yd are unconstrained by T along the y u = yd line, and the limits from 5* are less 
severe, the only bound we can obtain on y comes from the most stringent bound on yi due to 
the T parameter. We then see that typical values M ~ 1, 1.5 TeV and y — 1 are permitted 
for n v = 2. For larger values of n v , the allowed split between the mass eigenstates is smaller, 
but still large enough that a coupling y ~ 0.5 is not unreasonable. 

To summarize, we conclude that for a sufficiently large range of the parameters the fit to 
electroweak observables is in agreement with the existence of extra vector-like fermions. This 
parameter space is characterized by 

p r wl: y< 1.2 -1.6 

Pd = Pu = p q > Pi =Pe ■ y< 0.6 -1.2 (4.21) 
together with the current direct search limits 

M u M- e > 200 GeV 

M q , M a , M- d > 250 GeV. (4.22) 

5 Numerical analysis 

An interesting feature of the pure-glue hidden valley is that the v-glueballs can have many 
decay modes, depending on their quantum numbers and on the values of the various parame- 
ters. To give a general survey of decays is beyond the scope of this work, since a more detailed 
treatment would have to incorporate precise values of the matrix elements which at present 
are unknown. However, a few examples are useful to illustrate the main qualitative features 
of v-glueball decays. 

5.1 Decay patterns of v-glueballs 

The relevant parameter space consists of the mass scale M, the ++ mass mo, the coupling y 
and the Higgs mass rrtu- Here we present our results in terms of tuq which is more transparent 
than the confining scale A^, since mo is the relevant parameter for LHC studies. To simplify 
the discussion in the SM , it is convenient to assume that the Higgs boson is SM-like and fix its 
mass to be m# = 120 ("low mass" range) or m# = 200 GeV ("high mass" range). Moreover, 
since the branching ratios depend on M and y only through the combination yM, our problem 
is reduced to a two-dimensional parameter space described by yM and m . 
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For our estimates below, we will use the lattice results [T3] 

47ra„F^ ++ = 3.06mg, Fj ++ = 0.03m^ (n„/3), 47ra„F^_+ = 0.83mg. (5.1) 

We also need estimates of the other v-glueball decay constants and transition matrix elements, 
which at this point are unknown. A reasonable educated guess is that they are of order F^ ++ , 
which is the only known decay constant that is largely independent of the size of the v- 
glueballs. We can also guess M"__ 0++ ~ 1/Fq ++ , as it would be true for pion emission. In 
the following, we will therefore assume 

= ^M 2 s « o++ m * FJ ++ , M^ 1+ _m 1+ - = ^M?__ 0++ ~ ^m G /F s 0++ . (5.2) 

Later in this section, we will comment on how large deviations from this guess might affect 
our estimates. 

Using the decay rates expressed in section 3, we identify r^ 6 ^ and as the summed contri- 
butions to the decay rates from dimension-six and dimension-eight operators, respectively. The 
corresponding branching ratios are denoted by BR^ and BR^ 8 \ satisfying BR^ +BR^ = 1. 
To illustrate the dependence of the branching fractions on yM and m , we present in figures H] 
and [5] contours of constant BR^ in the uiq vs yM plane for various choices of v-glueball states 
and mu = 120 GeV, 200 GeV. In figure HI we see that the values of BR^} + are large over 
most of the parameter space plane, except at very low yM. By contrast, the values of BR,f} + 
are small, except for a small region at large yM. The most interesting decay pattern is found 
in the 1 v-glueball. We see that BR^}_ typically dominates when mi — — m\+- < m# and 

the emitted Higgs boson is off-shell, whereas BR\1_ dominates when m\ — — m\+- > m# and 
the emitted Higgs boson is on-shell. An analogous behaviour is found for inn = 200 GeV, as 
demonstrated in figure [5j 

The values of yM determine to a large extent the decay pattern of v-glueballs with three 
distinctive energy regions: (i) yM <C 1 TeV, (ii) yM 3> 1 TeV and (Hi) yM m 1 TeV. Recall 
that these ranges of parameters are further constrained by existing experimental limits, as 
described in section 4. As will be discussed below, the common feature of all these differ- 
ent regimes is the great diversity of decay modes and lifetimes, even for a given choice of 
parameters. 

(i) yM 1 TeV. For sufficiently small yM, the contribution of the dimension-six opera- 
tors can be effectively disregarded and, hence, the dimension-eight operators have the largest 
effects on v-glueball decays. As shown in [UJ, annihilation decays dominate the ++ , 2 ++ , 

8 Here the coupling constant Aira v is included alongside Fq ++ and Fq_ + to make them renormalization 
invariant so that there is no question at which point Aira v is normalized. On the contrary, Fj ++ is renor- 
malization invariant as is, since it is the matrix element of the energy-momentum tensor which is known to 
be scale invariant. Also, since the values reported in [13j are not expressed in a continuum renormalization 
scheme, we have converted <? 2 Fj ++ to Fj ++ using the value of the lattice parameter j3 = 6/g 2 = 3.2 as an 
approximation to the continuum limit. 



24 



200 400 600 800 1000 

m (GcV) 



200 400 600 800 1000 

TOo(GeV) 



200 400 600 800 1000 

roo(GeV) 



Figure 4: Curves of constant branching ratio BR^ in the parameter space (mo,yM) for 
various representative states and m# = 120 GeV. Left: ++ , Center: 2 ++ , Right: 1 . 
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Figure 5: Same as figure H] for m# = 200 GeV. 



0~ + and 2 h v-glueballs. Their branching ratios are dominated by decays to gg, with decays 
to 77 having a branching fraction of ~ 0.4%, assuming the X fields form complete SU(5) 
multiplets of equal mass. If the colored particles are much heavier than the uncolored parti- 
cles, then decays to electroweak bosons and photons can dominate. Most other states decay 
by radiatively emitting a photon, or a Z boson. In addition, the 1 state can also decay to 
standard model fermions via an off-shell 7 or Z. The 3 ++ state is special in that its three- 
body decay mode 3 ++ — > ++ gg could be of the same rate as the radiative 3 ++ — > l +_ 7 decay. 
The clearest signatures for this regime are likely to be the two-photon resonances from the 
prompt annihilation decays of C-even v-glueballs. In addition, for part of the parameter space, 
the v-glueballs are rather long-lived particles, so that they produced displaced vertices in the 
detectors, which can be an experimentally challenging, but certainly important signature. 

(ii) yM 3> 1 TeV. In this case, the Higgs couplings to v-glueballs dominate over all other 
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couplings. As a result, the decay pattern is relatively simple. In the PC = ++ sector, we 
find that the ++ v-glueball annihilates directly into pairs of standard model particles via 
++ — > h* , with the same final states and branching fractions of a Higgs boson with mass 
mo- For > 2rrih, the mode ++ — > hh opens up, in addition to ++ — > h*, contributing 
20 — 25% to the total width with the decays into weak boson and top quark pairs accounting 
for the remaining 75 — 80%. The 2 ++ v-glueball decays predominantly to the lighter ++ 
by radiatively emitting a Higgs boson. In the PC = — h sector the lightest states are the 
0~ + and 2~ + v-glueballs. The 2~ + v-glueball decays in a similar way to the 2 ++ , with the 
rate dominated by 2~ + — > v h. Although the 2~ + v-glueball is heavier than the ++ , the 
2~ + — > ++ h decay is not allowed, due to parity and conserved angular momentum. The h 
v-glueball is a special case. Without explicit breaking of parity, the _+ v-glueball decays 
slowly into gg via dimension-eight operators. 

Turning to the PC — H — sector, we find that most of the states can decay by emission of a 
Higgs boson, through processes of the form K — > Q K >h, where 9 K , Q K > denote two v-glueballs 
with given quantum numbers. Since C-changing transitions are not allowed by the 5* operator, 
the v-glueballs in that sector typically undergo a cascade decay, radiating a Higgs boson at 

each step, which ends at the lightest l + ~ state. The v-glueballs in the PC = sector 

decay in a similar manner, with the lightest 1 v-glueball decaying to the l + ~ v-glueball 
with the emission of a Higgs boson. Since the l +_ v-glueball is not permitted to decay via 
dimension-six operators, its dominant decay modes are photon-radiative transitions to C-even 
v-glueballs, as in (i). 

We should note that this theoretical regime may be difficult to access in practice. For 
the 2 ++ v-glueball (and hence the 2 _+ ), this requires yM > 10 TeV (see figures 0] and [5]) . 
As shown in f)4.2ip . the heavy mediators can at most have Yukawa couplings y ~ 1 in order 
to avoid potentially dangerous electroweak corrections. With M ~ O(10 TeV), the rate for 
production of v-glueballs (oc a gg ^xx *C lfb) is rather small, rendering this region of parameter 
space experimentally inaccessible to the future experiments at the LHC. 

(Hi) yM ?a 1 TeV. For Higgs couplings with intermediate strength, the v-glueball decay 
pattern is more complicated because the couplings to the gauge bosons become important, 
leading to a interesting interplay between dimension-six and dimension-eight operators. From 
the point of view of prospective experiments at the LHC, this is also the most interesting 
regime, due to the diversity of v-glueball decay channels and variability of lifetimes as well. 

The ++ v-glueball still decays predominantly to pairs of standard model particles via 
++ —> h* or, if kinematically allowed, ++ — > hh. The branching ratios for the main decay 
channels of the ++ v-glueball are shown in figure [6] for yM = 1 TeV in the cases of run = 
120 GeV and m H = 200 GeV. 

For the 2 ++ and 2 h v-glueballs, the dominant decay mode is gg. This is because the 2 ±+ — > 
±+ h(h*) decay is phase-space suppressed by a 10~ 8 — 10~ 3 factor. For m 2 ++ — m ++ > ran, 
the phase space suppression is smaller ~ 0.01, but in this regime the (m2++/M) 4 suppression 
of D = 8 operators is inefficient, rendering the D = 6 operators subdominant. 
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Figure 6: The branching ratios of the ++ v-glueball as a function of mo for yM ~ 1 TeV. 
Left Panel: ran = 120 GeV. Right Panel: run = 200 GeV. For clarity, only the main decay 
modes are shown. 



The dominant decay modes of the 1 v-glueball are transitions with a Higgs boson in the 
final state, 1 — > 1 H h, for m\ — — m^- > itlh or, photon-radiative decays to the lighter 
v-glueballs in the C-odd sector for m 1 — — m^- < The decays of the other v-glueballs 
in the C-odd sector proceed in a similar fashion. In addition, the decay of the 1 v-glueball 
into fermion pairs can play a significant role for mi — — mi+- < m#. 

Since the l + ~ and 0~ + v-glueballs cannot decay via their couplings to the Higgs boson, 
their dominant decay modes are gg for the _+ v-glueball and radiative transitions to the 
C-even states with emission of a photon for the l +_ v-glueball. 

The 3 ++ state is more complicated. With many contributing decay channels and unknown 
form factors, it seems impossible to estimate which decay mode is dominant. Indeed, a simple 
estimate suggests that the 3 ++ — > ++ gg, 3 ++ — > l + ~7 and 3 ++ — > 0~ + h decays are all of the 
same order. 

Summarizing all three cases, we expect the v-glueball decays to produce any of the kine- 
matically allowed final states of the Higgs boson, such as bb, t + t~, W + W~, etc, as well as 
possible photons, both singly and in pairs, and gluon pairs. If kinematically allowed, multiple 
production of standard model Higgs bosons from cascade decays can also proceeds with a 
sizable rate. A typical final visible state would then be of the form bbbb, WWWWbb, bbr + T~ , 
bbgg, 6677, and so on. 



5.2 Lifetimes 

Meanwhile, with so many v-glueball states and decay channels, the lifetimes of the v-glueballs 
can vary over many orders of magnitude. This is already clear from the fact that the lifetimes 
are very sensitive to both M and itiq. In our approach, the large contributions of the dimension- 
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six operators also suggest a wider spread of the lifetimes than that found when the Higgs 
couplings are absent, as in [TT]. In figure [TJ we plot the lifetimes of some of the v-glueballs as 
a function of m for the three representative regimes studied in the last subsection. We see 
that for yM ~ 1 — 10 TeV the various v-glueball states have lifetimes that typically span 5-6 
orders of magnitude, for any given choice of the parameters. 

An important consequence of the large spread in the lifetimes is that there is a significant 
probability that one or more of the v-glueballs will often decay a macroscopic distance away 
from the primary interaction vertex. Only one of these states needs to be both long-lived and 
frequently produced to provide a strong signature of new physics. 

If production rates are substantial, displaced vertices are expected for average decay 
length^ ct of the order 10 -4 — 100 m, resulting in lifetimes of 10~ 12 — 10~ 6 sec. Outside 
this range, the displaced-vertex signature is no longer present, either because the lifetime is 
much longer than 10~ 6 sec, so that a typical v-glueball produced at the LHC will escape 
the detectors, or because the decays are prompt and the v-glueball decay lengths cannot be 
resolved. 

As shown in figure [7J displaced vertices could be observed for 10 GeV < m < 200 GeV 
and < yM < 1 TeV or for 10 GeV < m < 400 GeV and yM > 1 TeV. As long as m 
becomes greater than about 400 GeV, all v-glueball decays will be prompt and the displaced- 
vertex signature will be absent. In the case yM « the l +_ and 1 v-glueball lifetimes reach 
values high enough for displaced vertices for mo between 50 GeV and 150 GeV. The ++ , 0~ + , 
2 ++ and 2~ + v-glueballs may also decay with displaced vertices in the very low mass range 
m o ^ 70 GeV. On the other hand, in the case yM ^10 TeV the spread in the lifetimes is two 
or three orders of magnitude larger, so even larger v-glueball masses 150 GeV < mo < 400 GeV 
may allow for displaced vertices. This is the case of the 2 ++ and 2 _+ v-glueballs, as shown 
in the middle panel of figure [7J Their lifetimes are so long that these v-glueballs may escape 
the detectors for m < 100 GeV. In the low mass range 50 GeV < m < 200 GeV other 
v-glueballs such as the l +_ and 1 also have a significant chance to decay with displaces 
vertices. In this case, the ++ and 2 ++ v-glueballs will remain short-lived over most of the 
parameter space down to masses of order 50 GeV where displaced vertices may start to occur 
for these states as well. Finally, the case yM « 1 TeV shown in the right panel of figure [7J is 
an example of an intermediate situation, with the ++ v-glueball being short-lived, except for 
very-low masses m < 50 GeV, and at least five states having a chance to decay with displaced 
vertices in some part of the parameter space: the 1 H and 1 v-glueballs for m < 200 GeV, 
and the 2 ++ , 2"+ and 0" + v-glueballs for 20 GeV < m < 100 GeV. 

From the analysis above it is quite evident that long-lived resonances are a common feature 
of pure-glue hidden valleys for an ample range of parameters. Nonetheless, detecting long- 
lived particles presents several experimental challenges. Displaced jets, though they have 

9 Note that the actual v-glueball decay length I in the laboratory frame, distributed according to P(l) oc 
e -i/icr w j^ n t}j e Lorentz boost factor 7 = 1/ \J\ — fi 2 , may be actually much smaller than the average decay 
length ct. So even if cr far exceeds the dimensions of the detector, there is still a significant chance that these 
v-glueballs will be observable, providing the rates are substantial. 
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Figure 7: Lifetimes of the v-glueballs as a function of the v-glueball mass scale mo for three 
representative regimes: yM ~ (Left panel), yM 10 TeV (Middle panel), and yM f=s 1 TeV 
(Right panel). 

no standard model background, suffer from substantial detector backgrounds from secondary 
vertices inside the detector material. On the other hand, displaced leptons are technically 
easier, and have much less background, but branching ratios to leptons are usually quite 
small. Displaced photons also present several experimental challenges, since one has to detect 
not only the position where the photons enter the electromagnetic calorimeter but also their 
angle of incidence. Thus the issue of detecting displaced vertices is not straightforward, and 
requires novel analysis strategies. 

Interestingly enough, the v-glueball masses can be as low as a few GeV's and at the same 
time the lifetimes could be short enough that a few v-glueball events could be observed at the 
LHC Specifically, let us examine closely the case of the scalar ++ v-glueball. From ( 13. 2 p we 
see that for X particles of order 1 TeV the lifetime of the ++ v-glueball is of order 

\ mo J \ m„ J \ ml) 

We see that for m > 10 GeV the ++ state can decay inside the detector. It is important 
to remark that even for cr well above 10 meters, a significant fraction of the v-glueballs will 
still decay inside the detector. If detector backgrounds are sufficiently low, and production 
cross-section is substantial, discovery would also be possible, in principle, for states in the 
low-mass range 1 GeV < mo < 10 GeV. 



5.3 Uncolored X particles 

If there is a substantial hierarchy between the colored X particles and the uncolored X par- 
ticles, the latter may dominate v-glueball decays. Decay rates to a pair of SM gluons in the 
C-even sector are proportional to the coefficients xi (EH2J), which would be suppressed. If 
the particles carrying QCD color are heavier by factors of 2 or more, the hierarchy of decay 
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rates changes, as X3/X7 ~ 10 2 — 10 4 . In this case, decays to gluons no longer dominate 
the partial widths for the C-even states. Also, the regime for m and yM in which 
dominates over r^ 6 ^ is somewhat reduced. Without any detailed computation, we may infer 
from figures H] and that the decays of the ++ v-glueball, with its large branching fraction to 
off-shell Higgs boson, will not be much affected. For tuq- < 2mw, the dominant decay mode 
of the _+ v-glueball is now 77 with a branching ratio of 90% with the decays into gg and Z7 
accounting for the remaining 10%. In the high mass range m - > 2mw, the 0~ + v-glueball 
decays dominantly into WW with a branching ratio of ~ 70% followed by the decays into 
ZZ, 77 and Z'j with branching ratios of 20%, 5% and 5%, respectively. Also, the lifetime of 
the _+ v-glueball is longer. For the 2 ++ and 2~ + v-glueballs, the decays into photon pairs 
typically dominate for m 2 ±+ — m ±+ < mn, while the Higgs-radiative transitions 2 ±+ — > ±+ h 
can dominate for m 2 ±+ — m ±+ > ran- The decays of the v-glueballs in the C-odd sector are 
largely unaffected. We will comment further in the conclusions and in our LHC study [3D] , 
but suffice it to say that the fact that the annihilation into photon pairs may dominate the 
lifetimes of some v-glueballs (the h , 2 ++ , 2~ + ) has important consequences for experimental 
searches. 



5.4 Uncertainties in the v-glueball matrix elements 

In another respect, we should emphasize that even though our estimates are in general robust, 
they are subject to significant uncertainties, due especially to the many unknown v-glueball 
matrix elements. Although lattice computations are at present not available, it is of interest to 
see what one might expect once lattice computations are incorporated. Of course, the decays 
of the ++ v-glueball will remain unchanged, since its branching ratios are independent of 
the decay constant F^ + . The same is true for the l +_ and h v-glueballs. On the contrary, 
the decays of the 2 ++ and 1 v-glueballs depend on ratios of the unknown non-perturbative 
matrix elements, 

M 2++o++ m o Mf__ 1+ _m 1+ - 

r 2 + = — r i- = ~^rm ■ ( 5 - 4 ) 

r 2++ lvl l — 0++ 

Since there is no reason for the form factors M2 ++0++ and 1+ _ to be significantly sup- 
pressed or enhanced at k 2 = m 2 H , the ratios r 2 + and r x - are probably of order 0(1) (actu- 
ally somewhat smallei0, especially at large n v ). To illustrate the possible range of impact 
lattice computations may lead to, we can simply make the rescaling yM — > t/My^y and 
yM — > yM^Jri- in figure SJ One can see that when r 2 + is varied by a factor of 10 in both 

directions around r 2 + ~ 1 for yM ~ 1 TeV, the ratio BR^l changes by around 10 — 20%, and 
the dominant mode of the 2 ++ v-glueball is still gg. A similar conclusion is obtained for the 
2~ + v-glueball. For the 1 v-glueball, in contrast to the 2 ++ and 2~ + , the branching ratios 
are more uncertain. The reason is that the contributions from D = 6 and D = 8 operators are 
comparable in strength, so a small enhancement in one of the matrix elements can lead to a 



10 Using standard large-n„ counting rules and taking into account explicit factors of in the coupling 

constant, one obtains r 2 + ~ 0(l/n v ) and r±- ~ 0(l/y/n^); see 
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large effect on the branching ratios. For T\- ~ 0.1, the photon-radiative decay gets enhanced 
very greatly and could dominate the branching ratio in most of the mass range, while r\- ~ 10 
tends to enhance the Higgs-radiative decay. Therefore, additional lattice computations would 
be needed to distinguish between these two modes. Fortunately, it turns out that the impact 
of precise lattice computations on the phenomenology of the lightest v-glueballs ++ , 2 ++ , 
0~ + and 2~ + , which, as will be argued in our LHC study, are probably the most copiously 
produced, is in general very mild. 



6 Other extensions 

In this section, we consider two simple extensions that may alter the phenomenology of v- 
glueballs in a number of ways, with a special emphasis on the h and l + ~ v-glueballs. 
Firstly, we analyze the two Higgs doublet model (2HDM), which is the simplest extension of 
the standard model Higgs sector, and secondly, we study the possibility of CP violation in 
the theory and its implications for v-glueball decays. We will not attempt to systematically 
compute all the decay rates, but simply point out a few salient features of these extensions. 



6.1 2HDM 

As in the MSSM, we consider the SM with two Higgs doublets H u and H d , where H u only 
couples to up- type quarks and neutrinos and Hd only couples to down-type quarks and leptons. 
In this model, there are five physical Higgs bosons: two charged scalars (iP 1 ); two neutral 
scalars (H and h); and a neutral CP-odd scalar (A). The presence of the CP-odd scalar A 
is of particular interest. Since the A is CP-odd, a gauge invariant CP-conserving coupling 
to the v-gluons must be of the form e^ ual3 tr T^T^A^, with <3> = l,H,h. Because of this 
coupling, new possible decays of the pseudoscalar 0~ + state emerge in the 2HDM. These are 

f, (6.1) 

(6.2) 
(6.3) 
(6.4) 

For all other v-glueballs, there are new contributions to their decays in the 2HDM, but the 
decay pattern as described in sections 3 and 5 remains basically unchanged. For the ++ 
v-glueball, the new channels include ++ — > hh, Hh, AA as well as the off-shell Higgs decay 
++ — > h*. For all other states, new decay modes arise from the processes 6 K — > Q K >4>, where 
now <p = H,h, A. Recall that in the single-Higgs model, the lowest-dimension operators that 
can induce decays of the l +_ v-glueball arise at mass dimension eight. A similar conclusion 
applies to the l +_ v-glueball in the context of the 2HDM. Therefore, while the CP-odd scalar 
A makes a drastic change on the 0~ + , it does not affect the l + ~ v-glueball. 
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To estimate the effect of the CP-odd scalar A on the branching ratios, we will need to find 
the effective Lagrangian. Using the basic results of section 2, the effective Higgs couplings to 
v-gluons, induced at one-loop by the X particles, is given by 



£(6) 



Ot v 
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vlHiH u + ylH\H d 



tr J~ uvJ~ 
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H u — HI 
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Here y u and yd represent the couplings of the X particles to the Higgs doublets H u and H d , 
respectively. The couplings of the physical Higgs bosons to v-gluon pairs are readily obtained 
by including the appropriate mixing angle factors 
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where the coeffients ^ and & are given by 

£hh = y 2 u s 2 a + y 2 d c 2 a ihh = y 2 u c 2 a + y 2 d s 2 a £ Hh = 2s a c a {y 2 u - y 2 d ) 

Uh = yls a cp - y 2 d c a sp £ Ah = y 2 u c a cp + y 2 s a sp £ A a = yl<?p - y 2 d s 2 p 

iH = y 2 u s a sp + y 2 d c a c p £ >h = ylc a sp-y 2 d s ol c j} £ A = spcp{y 2 u - y 2 d ) 

where a the mixing angle and the ratio of the vacuum expectation values tan (3 = v u /v d . The 
combination v\ + v 2 is fixed by the electroweak scale vh = \f^ d + v 2 = 246 GeV. 

The decay rates can be readily extracted from our general results in section 3 by replacing 
y 2 with the appropriate couplings in (16.71) . For the scalar ++ state, the decay rates for the 
off-shell Higgs decay ++ — > H*/h* — > can be read from (13.21) through the substitution 
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while the rate for the two-body decays ++ — > HH, hh, Hh, AA is obtained from (13. 3p by 
replacing 
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where i,j = HH, Hh,hh, AA and are the cubic self-couplings in the 2HDM. Likewise, 
the width of the decays 6 K — > Q' K &i, = H, h, A, can be obtained from (13 .9p and ( 13. 1Q|) by 
substituting y 2 — > £j. Finally, we summarize the rates for the new decay modes of the h 
v-glueball fl6.lU6.4j) . 



> ff ~ {2nM 2 (m 2 A - w 2 ) > h ^f [m ° 



(6.10) 
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2 2 2 \1 1/2 

m H ,m A ;m -)J 

//■ "4; m o-)] 



3m 2 z £ A 



2« 



[g(m 2 h ,m 2 A ;m 2 - 



)] 



1/2 



(6.11) 
(6.12) 
(6.13) 



The expressions above are clearly very model dependent. First, notice that for y u ~ y^, £4 
is very small, and the decay rate of the CT + v-glueball is suppressed. On the other hand, for 
y d ~ 0, a quick check shows that this rate, whose ratio to the ++ width is (for small masses 
m < m h ) 

r -n 



r ++ 



0.5 



m\ m 

S rv ~ ~\~ Cq 



'a o 

mi 



"H m h. 

is not negligible. This ratio may be anywhere from ~ 0.1 up to about 4 depending on the 
mixing angle and the masses of the scalars. Unless the h decay rate is unduly suppressed, 
the lifetimes of the ++ and 0~ + v-glueballs are in general within one order of magnitude from 
each other. 



(6.14) 



6.2 Models of explicit CP violation 

It is of interest to see what one might expect once CP- violation effects are incorporated in 
the theory, in the context of a single-Higgs model. For example, parity-violating interaction 
terms can be induced in the effective Lagrangian if we assume that the couplings of the heavy 
particles X r to the Higgs boson are complex. In this case, a P-odd C-even interaction is 
generated of the form if'iftr T llv T iXV . As in the 2HDM, this interaction drastically changes 
0~ + decays, since without it, the 0~ + v-glueball would decay via dimension-eight operators, 
suffering an extra suppression. The 0~ + v-glueball can then decay via h — > h*, producing 
any of the kinematically-allowed final states of the Higgs boson. 

Similarly, one can contemplate interactions that explicitly break C-invariance. Interestingly 
enough, the lowest-dimension operators of this kind (e.g. H* Hd a i >e f c deJ~% l/ J~, b M ^~aB^'aB) ar i se 
at dimension ten and, therefore, their effects are extremely suppressed by extra powers of 
1/M. The dominant decay modes of the l + ~ v-glueball are then induced by dimension-eight 
operators, even if we allow for CP-violation. 

In the rest of this section, then, we will focus on P- violating but C-conserving interactions. 
To estimate the decay widths, let us consider the following Lagrangian, in two-component 
notation, 

Cmass = VdX q H^X c d + y d X c q HX d + h.c. (6. 15) 
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where yd is a complex Yukawa coupling^ By integrating out the X q and X d particles, we obtain 
the following dimension-six operator 



11. 



where y = Re yd and y = Irayd- Here T^y = (l/2)e flu x a J-'' 1L/ denotes the dual of the field 
strength tensor. The first term on the right-hand side of (16.161) is the same operator that we 
had already found in (12. 2p . The second term of (I6.16P is new and violates P. As mentioned 
above, the P-odd interaction allows the pseudoscalar 0~ + state to decay into SM particles via 
s-channel Higgs-boson exchange 0~ + — > h* — > ((, where ( denotes a standard model particle. 
The width of the decay is given by 

f 2yyv H a v F^ \ 2 SM , 
\irM z {m H — m^J / ss 

where F^_ = a v (0|tr J r /ii/ .F At!/ |0 ++ ) is the 0~ + decay constant. The same operator also induces 
the decay CT + — > ++ h, with partial width 

r °- + - 0++/l " 16^7 ViXP ) [9( m H,m 0+ ,ra -)\ (6.18) 

where now Mg_ 0+ is the transition matrix. This decay is phase-space suppressed for m -+ — 
m ++ < m h , but may be quite significant for m -+ — m ++ > m h when the radiated Higgs 
boson is onshell. For example, for ~ 100 GeV we obtain BR(0~ + — > ++ h) ~ 0.3. 

In the SM, the most stringent limits on the size of the CP violation effects come from 
experimental limits on the electric dipole moment of the neutron. The same experimental 
constraints also allow us to place limits on the CP- violating operator in (I6.16p . To see this, 
let us recall that the X particles can also carry QCD color. Then, the following dimension-six 
operators are induced: 

£ (6) = a " { f~J 2) H^H tr G m GT + ^rj^H^Htr G^. (6.19) 

where now G^ v denotes the SM gluon field strength tensor and G^ v its dual. The CP violating 
term in (I6.19P with the Higgs field replaced by its vacuum expectation value contributes to the 



11 The effective Lagrangian (|6.16l) can be evaluated either by following the procedure displayed in appendix 
A or, equivalently, in the limit of vanishing Higgs momentum by taking derivatives of the v-gluon self-energy 
and the axial anomaly: 

= H^H/v 2 H (yjd 2 /dy 2 d + yfd 2 /dyf) C, 

where 

C = (a v /6n) \n{k 2 uv /AetMHl) tr T^T^ + {a v /2m) \n(detM/detM*) tr T» v 



with M 



M y d v H 
y d v H M 
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9 angle. We assume that this contribution to the 9 angle is removed by whatever mechanism 
that solves the strong CP problem in QCD. The next operators in the expansion around the 
Higgs VeV contain one or two powers of the Higgs field. Their contribution to the electric 
dipole moment of the neutron can be estimated using Naive Dimensional Analysis [IT] W2\ . 
Following the analysis of [13| 01] , we estimate the neutron electric dipole moment as 



*■ ~ W ~ ^ — " ^ « <* TeV/M 2 ). (6,0) 

where ln]^ ~ 1190 MeV is the chiral-symmetry-breaking scale. Combining this result with 
the current experimental bound on d n , we derive a limit on 2yy 

\2yy\ < 10 x (M/l TeV) 2 . (6.21) 

Consequently, values for \y\, \y\ in the range 0.01 — 1 are consistent with the existing experi- 
mental limits. 



7 Conclusions 

In this work, we have investigated a particularly challenging hidden valley scenario with a 
broader class of couplings to SM particles than those considered in [JTJ. In particular, we 
have focused on the effect of dimension-six operators by which the hidden sector interacts with 
the standard model through the Higgs sector. 

The resulting v-glueball phenomenology is fairly complex in this scenario, but there are 
some simple features. In particular, we find the following interesting signatures: 

• Decays of the ++ v-glueball through h*, producing any of the kinematically-allowed 
Higgs final states such as bb, t + t~ , WW, etc. 

• Multiple Higgs boson emission, from cascade decays 6 K — > K >h, or annihilations ++ — > 
hh. 

• Due to the diversity of v-glueball states, and the presence of operators of different mass 
dimension, the lifetimes can vary at least over 5 or 6 orders of magnitude, for any given 
choice of parameters (see figure [7]). This has two immediate consequences: (i) Displaced 
vertices are quite common, which could potentially serve as a discovery channel, and (ii) 
v-glueballs can be as light as a few GeV's and still be visible. 

• For sufficiently small v-glueball masses, a different opportunity arises in the form of 
non-standard Higgs decays such as h — > with branching ratio to v-glueballs of 
order 10" 3 - 10~ 2 . 

In addition, other possible final states include gg, 77, and radiative transitions with emis- 
sion of a photon, all of which are induced by dimension-eight operators [IT] . The rare diphoton 



35 



channel is particularly interesting because, thanks to its moderate QCD background, it may 
serve as the discovery channel for the v-glueballs at the LHC. 

We should nevertheless emphasize that some of our results are subject to significant the- 
oretical and numerical uncertainties. In this respect, it would be interesting to know the 
spectrum of pure-Yang-Mills theory for gauge groups other than SU(3), as well as the various 
v-glueball matrix elements that arise. With enough motivation, such as a hint for discovery, 
these could in principle be determined by additional lattice computations. If the v-sector gauge 
group is not SU(n v ), some of the v-glueballs may not be present; for instance, for SO(n v ) or 
Sp(n v ) gauge groups the C-odd sector is absent or heavy. However, as explained in [11], the 
lightest C-even v-glueballs are expected to be present in any pure-gauge theory, with similar 
production and decay channels, so at least for them, the basic features of a pure-gauge hidden 
valley are expected to be retained. We also have not considered higher-order corrections to 
the decay rates, and they should be taken into account when precise predictions are required. 

We have seen that the operators considered in this paper are not heavily constrained by 
current experimental searches or precision electroweak data, thus leaving ample parameter 
space to be explored by the next generation of hadron collider experiments. Application of 
our results for phenomenological studies, particularly as relevant for the LHC, will be carried 
out in a companion paper [10] ■ We expect that detection should be feasible, if the mass M of 
the X particles is small enough to give a reasonable large cross-section, and A„ is large enough 
to ensure the v-glueballs decay inside the detectors. Decays to bottom quarks and gluons 
often dominate, but they suffer from a huge multijet background. At the Tevatron and at the 
LHC, the v-glueballs are most likely to be found in searches for diphoton resonances in events 
with 2 photons plus jets or with 3 or more photons and possibly additional jets, or in searches 
for displaced decays to jet pairs, W/Z pairs, or photon pairs. A potentially novel signature, 
which could be currently searched for at the Tevatron, arises in the form of two b-tagged jets 
plus diphoton events. Contrary to the models of [15J HOI H7], the branching ratio to 6677 
in our model can be 0(1), since the bb and 77 final states originate in different resonances, 
leading to a potentially discoverable signal despite the background of 2j + 27. Altogether, 
these signatures provide an interesting opportunity to search for new physics at the Tevatron 
and at the LHC, which motivates further theoretical and experimental study. 
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A Computation of cfL 

Definitions 

In the following the momenta of the initial gluons (pi,P2) and final Higgs bosons (ki, k 2 ) are 
all chosen to be incoming. We make use of the following invariants 

s = (pi + P2) 2 t = ( Pl + k 4 ) 2 u=( Pl + k 3 f (A-l) 

obeying the condition 

s + t + u = kl + kl (A-2) 
where we have assumed that the gluons are massless, p\ = p\ = 0. 

Pasarino-Veltman functions Co and Do 

The computation of the Feynman graphs is performed using the Passarino-Veltman decom- 
position of the one-loop tensor integrals [48]. For our results below we will make use of the 3 
and 4 point scalar functions Co and D$. The explicit expressions are rather lengthy and have 
to be handled numerically. However, one can derive an integral representation 



i>1 px 

C (n,r 2 ,s,M 1 ,M 2 ,M 3 ) = / [~r 2 x 2 - ny 2 + (s + r x + r 2 )xy+ 

Jo Jo 

+ (M 3 - M 2 + r 2 )x + {M 2 -M 1 + s- r 2 )y - Ms]" 1 dydx (A-3) 



pi px py 

D (r 1 ,r 2 ,r 3 ,r 4 ,s,t,M 1 ,M 2 ,M 3 ,M i )= / / [-r 3 x 2 - r 2 y 2 - r lZ 2 + 

Jo Jo Jo 

+(—t + n + r A )xy + (s + t — r 2 — r 4 )xz + (— s + r x + r 2 )yz + (M 4 - M 3 + r 3 )x+ 



+ (M 3 - M 2 + t - r 3 )y + (M 2 - M 1 + r 4 - t)z - M 4 ]~ 2 dz dy dx (A-4) 



Fermion loop 



The amplitude for the fusion of two massless v-gluons into a pair of Higgs boson fields arises 
at lowest order in perturbation theory from a loop of heavy particles. The relevant graphs are 
depicted in figure [HJ We assume that there is a pair of vector-like fermion fields, a doublet 
and a singlet of SU(2)i, transforming under the fundamental representation of SU(n v ). The 
doublet and singlet have masses M and m, respectively. Gauge invariance constrains the 
amplitude to be 

2 2 

Mo = |Jf (p£p£ - Pi ■ V29T) Fu + • ■ ■ (A-5) 
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Figure 8: Graphs contributing to gg — > HH 

where e* e 2 are the polarization vectors of the incoming v-gluons. The form factor F u is given 
by: 

F D = f(s, t, u, k%, kl, m, M) + f(s, t, u, fc|, kj, m, M) + (M <H> m) (A-6) 

where 

f(s,t,u,kl,k\,m, M) = 
- ^(k 2 3 + k\ - 2(M + m) 2 ) [(kl - u)C (0, kl u, M 2 , M 2 , m 2 ) + (A; 2 - i)C (0, A; 2 , t, m 2 , m 2 , M 2 ) 
+ (A| - u)C o (0, u, kj, m 2 , m 2 , M 2 ) + (A; 2 - t)C o (0, f, A; 2 , u, M 2 , M 2 , m 2 )] 

777,^ 

+ 8— C (0, 0, s, m 2 , m 2 , m 2 ) 

77T, 2 iW 

- 2 — (-2(m + M) 2 + k\ + kl + —s) \D {0, 0, fcf , fcf, s, t, m 2 , m 2 , m 2 , M 2 ) 
s m 

+D (0, 0, kl kl s, u, m 2 , m 2 , m 2 , M 2 )] + 

— 2 (2m 4 s + 4m 3 Ms + m 2 (s{AM 2 - A; 2 ) + 2A;^A; 2 - A:|(2A:| + s)) + 2mM(2M 2 s - 2fc 2 (fcf - t) 
s 

+2k 2 t - s 2 - 2st - 2t 2 ) + (2M 2 - k 2 z - kl)(sM 2 + k\k\ - k 2 k 2 4 ) 
+s(2m 2 + 2M 2 -kl- k 2 4 )p 2 t ) D o (0, kj, 0, k 2 3 , t, u, m 2 , m 2 , M 2 , M 2 ) + 4 (A-7) 

with p t = 2(p 1 k 3 )(p 2 k 3 )/(p 1 p 2 ) - kl 
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Large fermion mass limit and effective Lagrangian 



The form factor can be evaluated by taking the limit M, m 3> s,t, u, k\, k\. At leading order 
in the mass splitting of the heavy fermions M ~ m, the form factor reads 



(A-8) 



The matrix element for g a gb — > HH can be obtained in perturbation theory from the following 
non-renormalizable interaction 



2 

,(6) _ V &v tt\tt± -p -ryv 



(A-9) 



B Classification of v-sector operators 

The classification of the operators Ov in terms of irreducible representations of the Lorentz 
group was carried out in [28]. (See also |TT] for a review.) In this appendix we summarize the 
operators which are relevant for our work. Here S, P, T, L and f^ 1 ' 2 ) label the decomposition 
of each operator Ov into irreducible representations of the Lorentz group. Their expressions 
are shown in tables and along with the J PC states that each operator can create acting on 
the vacuum. From now on, we denote the operators O^, where £ runs over different irreducible 



Operator C| 


jPC 


S = tr T^T^ 


0++ 


P = ti T^f^ 


o-+ 


T a p = tr T ^xJ 7 g X — \ g a pS 


2++, 0++ 


— tr J 7 ayj 7 a/3 — 2 (9^aT u /3 + g u pT^ a ~ 9fj,pT ua — g ua T^p) 


2++, 2-+ 


— 12 {9tMa9vj3 — 9^9va)S + e^ va pP 





Table 5: The dimension d = 4 operators, and the states that can be created by these operators [28] . 
We denote = \ e^pF^ . 



operators £ 



S,P,T,L,- 



Operator 0\ 



jPC 



- 

,(2) 



n 



ir ~F ~F ~F a P 

= tr F"F^Fp u 



Table 6: The important d 
shown l28l. 



6 operators. The states that can be created by these operators are 



39 



C The form factors My Jt 

In addition to the momenta p and q, we introduce a pair of tensors e Mli ... lflJ (p) and e Mli ... iM , (g) 
to represent the polarizations of the spin-J and spin-J' states respectively. It will be convenient 
to define "reduced" polarization tensors by contracting some of the indices with the vectors p 
and q, 

</,,-,<, M = W-,M»<f' 7 - K+1 • • • <t J W-,m„(<?) = W",^)^'-^ 1 • • •^' / '- 

(A- 10) 

The invariant form factors are constructed from the reduced polarization tensors and the 
vectors p and q. For the simplest J — > case, the only form factor is obtained when all the 
indices of the polarization tensor are fully contracted with q, 

e(p)=e Ml ,.., M >)<r...^. (A-ll) 

This gives the matrix element 

(0,q\S\J,p) = M s 0J e fllt ...^ J (p)q^---q^ (A-12) 

where now Mqj is the transition matrix which depends on the transferred momentum. For the 
more complex J — > J' — 1 case, there are two form factors parametrizing the matrix element, 

^ 2 ■■■ ( f J ^(<?)> w 2 -,m» <T • • • <f J Uq)p u - (A-i3) 

In the reduced notation, these are simply e M e M = and e e. In addition, we can contract the 
polarization tensors with a Levi-Civita tensor as follows, 

e^%(p) p p ^(g) = e^%, M2 ... ^ (p) q^ ■ ■ ■ q?' p p q a ~eM (A-14) 

Collecting all terms, we obtain the following matrix element 

(1, q\S\J,p) = Mfje^.^ip) q^--- <f> e"(g) + Mf^ 1)M2 .., M » q^ ■ ■ ■ <f> e„(q)p»+ 

M5^Ww(p) ^ • ■■^PpQMq) (A-15) 

We next turn to the more complicated matrix element for the J — >• J' = 2 transition. 
Firstly, we define the auxiliary tensors 

eM = W>-,^(p) ^ 3 --- 9 MJ , e,(p) = <„,,,.,,,,(,>) q^ ■ ■ ■ q^ (A-16) 

and 

e M (g) = e^(g)^. (A-17) 
From €fj, u (p), €(i(p), and e M (g) we can form the following three invariants 

e M G»)e M (g), (e„(p) <f) (e„(g) p") (A-18) 
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In addition, contracting with the Levi-Civita tensor gives two more invariants, 



e^ P %(p) P P qMl), V P q<r~e u \q). (A-19) 

Therefore, there are five form factors in the J — > J' = 2 case. The general rule to find all the 
form factors in the case J — > J' case is now clear. Without loss of generality we can assume 
J > J'. First we construct the reduced spin-J tensor with J' indices, 

e„i,-,/v(p) = <,,. ,, A I>)<1"' ■ ■■q ilJ (A-20) 
so that we have two tensors with J' indices. Then we form the invariants 

e fllr .. !fljl (p)e^-^'(q), ••• , e» e"( 9 ), e(p)e(g), (A-21) 

for a total of J'+l invariant form factors. Secondly, there are contractions with the Levi-Civita 
tensor as follows 

e^p,P P q a e^ 2 ,.. (p)^'"^' (q), ■ ■ ■ , e^jfq'e^p) ? (q) (A-22) 

giving J' additional invariants. Hence, the total number of Lorentz invariant form factors is 
2 J' + 1 for the transitions J — » J', J' < J. To see that this contemplate all the possibilities, 
we can apply angular momentum counting rules. By composing the spin J and spin J' states 
we obtain total spin | J — J'\, ■ ■ ■ , J + J'. For a scalar operator, these have to be combined 
with total angular momentum L = | J — J'\, ■ ■ ■ , J + J' in order to give total spin 0. So for 
J' < J we get 2 J 7 + 1 form factors associated with the different values of L. 



i 


1 


2 


3 


4 


5 




9ap 


g ap qpp<T 


qaPpqpPa 


g/3p taazw P Z (f 


(1/3 Pp taozw P Z q W 




9ap 


q a p P 


^apzw P q 







Table 7: The auxiliary tensors M.( J,l \ Here q a = q a / \[q^ and p a = p a / \fjP. 

To summarize, we show the matrix elements up to spin 3 that are needed to compute the 
decays of the v-glueballs in figure [TJ In terms of the auxiliary tensors shown in table [7J the 
matrix elements read, 

M<& = q 7 e^(p) e"{q) M%% (A-23) 
Mi = q, qp e a ^(p) ~e"(q) M™ (A-24) 
M 30 = q a qp q-ye a ^(p) (A-25) 

M$=e#(p)i"{q)M§% (A-26) 



41 



M$=b^(p)e>(q)M<g> (A-27) 

M w = q a q p e a \p) (A-28) 

M$=f(p)e>{q)M™ (A-29) 

M 10 = q a e» (A-30) 

where we denote q a = q a / For the sake of brevity, we have not shown the seven matrix 

elements that correspond to the 3 — > 3 transition. These can be obtained along the same line 
that led to 
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